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AN L,(L p )-THEORY FOR THE TIME FRACTIONAL EVOLUTION 
EQUATIONS WITH VARIABLE COEFFICIENTS 

ILDOO KIM, KYEONG-HUN KIM, AND SUNGBIN LIM 


Abstract. We introduce an L q (L p )-theory for the quasi-linear fractional equa¬ 
tions of the type 

d^u(t, x) = a 4 - 7 (t, x)u x i x j ( t , x) + f(t, x, u), t > 0, x E H d . 

Here, a E (0, 2), p, q > 1, and is the Caupto fractional derivative of order a. 
Uniqueness, existence, and L q (L p )-e stimates of solutions are obtained. The 
leading coefficients a XJ (t,x) are assumed to be piecewise continuous in t and 
uniformly continuous in x. In particular a^ (t, #) are allowed to be discontin¬ 
uous with respect to the time variable. Our approach is based on classical 
tools in PDE theories such as the Marcinkiewicz interpolation theorem, the 
C alder on-Zygmund theorem, and perturbation arguments. 


1. Introduction 

Fractional calculus has been used in numerous areas including mathematical 
modeling [24] [44] , control engineering mm, electromagnetism mm, polymer 
science b 02 , hydrology 01 SB, biophysics 01122, and even finance ESI 02- 
See also 03122123 US and references therein. The classical heat equation dtu = 
Ait describes the heat propagation in homogeneous mediums. The time-fractional 
diffusion equation d“it = Ait, a £ (0,1), can be used to model the anomalous 
diffusion exhibiting subdiffusive behavior, due to particle sticking and trapping 
phenomena (see [251 ES])- The fractional wave equation d t u = Ait, a £ (1,2) 
governs the propagation of mechanical diffusive waves in viscoelastic media (see 
1 2.311. The fractional differential equations have an another important issue in the 
probability theory related to non-Markovian diffusion processes with a memory (see 

E2E3). 

The main goal of this article is to present an L q (L p )~theory for the quasi-linear 
fractional evolution equation 

d?u(t,x) = a 13 (t,x)u x i x j{t, x ) + b l (t, x)u x i[t , x) + c(t , x)u(t, x) + f(t , x, it) (1.1) 

given for t > 0 and x £ R rf . Here a £ (0, 2),p, q > 1, and denotes the Caputo 
fractional derivative (see (12.31) 1. The indices i and j move from 1 to d , and the 
summation convention with respect to the repeated indices is assumed throughout 
the article. It is assumed that the leading coefficients a i3 (t,x) are piecewise con¬ 
tinuous in t and uniformly continuous in x, and the lower order coefficients b l and 
c are only bounded measurable functions. We prove that under a mild condition 
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on the nonlinear term f(t,x,u ) there exists a unique solution u to (11.111 and the 
L 9 (L p )-norms of the derivatives D^u, |/3| < 2, are controlled by the L 9 (L p )-norm 
of f(t,x, 0). 

We remark that there are a few other types of fractional derivatives such as 
Ricmann-Liouville, Marchaud, and Griinwald-Letnikov fractional derivatives. These 
three fractional derivatives coincide with the Caputo fractional derivative in our so¬ 
lution space M“’” 0 (T) (see [39l[3Tj for the proof). 

Here is a brief survey of closely related works. In m EH] an L q (Lp )-theory for 
the parabolic Volterra equations of the type 


d_ 

dt 


Lu + f 

\ J — ( 


k\ (t — s)u(s , x)ds 


Au + f(t,x,u), t£ R, a;eR d 


( 1 . 2 ) 


is obtained under the conditions k\ (t) > ct “ for small t, cq > 0, a £ (0,1), and 



aq p 


(1.3) 


The results of El EH also cover the case Co > 0, however it is obtained only for the 
case aX {t, x) = S u with the restrictions a £ (0,1) and (11.31) . If p = q, an L p -theory 
of type O) with the variable coefficients a IJ (t,:r) is presented in [|B] under the 
condition that a*- 7 are uniformly continuous in (t, x) and limi^i^oo a ZJ (f, x) exists. 
In [47] an Z^-theory is obtained for the divergence type equations with general 
measurable coefficients. Also an eigenfunction expansion method is introduced 
in EH t° obtain ^-estimates of solutions of divergence type equations with C 1 - 
coefficients. 

For other approaches to the equations with fractional time derivatives, we refer 
to m for the semigroup approach, to [51 1M] for ^([O,T],X)-type theory, where 
X is an appropriate Banach space, and to [7j for BUCi-p([0, T], X)-type estimate, 
where \\u\\ BUCi _ ? (io,t],x ) = sup tg(0jT] t 1_/3 ||u(t)||x- 

Our result substantially generalizes above mentioned results in the sense that we 
do not impose any algebraic conditions on a, p , and q. The conditions (11.31) and a £ 
(0, 1) are used in EJG2II, an d the restrictions p = q and a ^ — p’ 2 ^ 1 ^ 

are assumed in [46] . More importantly, in this article the condition on the leading 
coefficients a l;, (t,x) is considerably weakened. In particular, a lJ (t,x) depend on 
both t and x and can be discontinuous in t. Recall that if p > 2 then among above 
articles only [46] considers the coefficients depending also on t, but the condition 
p = q and the continuity of a l ° with respect to (t,x) are assumed in |46| . 

Another significance of this article is the method we use. The results of [8] f33ll46] 
are operator theoretic, j21] is based on iZ°°-functional calculus, and the method of 
HUGH works well only in the Hilbert-space framework. Our approach is purely ana¬ 
lytic and is based on the classical tools in PDE theories including the Marcinkiewicz 
interpolation theorem and the Calderon-Zygmund theorem. We obtain the mean 
oscillation (or BMO estimate) of solutions and then apply the Marcinkiewicz in¬ 
terpolation theorem to obtain L p -estimates of solutions. To go from L p -theory 
to Lq(Lp )-theory we show that the kernel appeared in the representation of solu¬ 
tions for equations with constant coefficients satisfies the conditions needed for the 
Calderon-Zygmund theorem. Perturbation and fixed point arguments are used to 
handle the variable coefficients and the nonlinear term respectively. 

The article is organized as follows. Some properties of the fractional derivatives 
and our main result, Theorem l2.91 are presented in Section 2. The representation of 
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solutions to a model equation and an ^-estimate of solutions are given in Section 
3, and BMO and an L g (L p )-estimate of solutions to a model equation are obtained 
in Section 4. The proof of Theorem 12.91 is given in Section 5, and sharp estimates 
of kernels related to the representation of solutions are obtained in Section 6. 

We finish the introduction with some notation used in this article. As usual N = 
{1,2 ,---}, R d stands for the Euclidean space of points x = (a; 1 , x d ), B r {x) := 
{y G R d : \x — y\ < r}, and B r := B r ( 0). For multi-indices 7 = (71,7^), 
7 i G {0,1, 2,...}, x G R d , and functions u(x) we set 

du 

“ x * = 3 ? = DJu = Df ■ ... ■ D?u, 

X 7 = (x 1 ) 71 (x 2 ) 72 • • • ( x d ) ld , | 7 | = 7i H-1-7 d- 

We also use D™ to denote a partial derivative of order m with respect to x. For an 
open set Q C R d by C£°(f!) we denote the set of infinitely differentiable functions 
with compact support in U. For a Banach space F and p > 1 by L p (U,F) we 
denote the set of F-valued Lebesgue-measurable functions u on 9 satisfying 


IMlL p (n,F) 



II p F dx 


1 Ip 


< 00. 


We write / € L Pt i oc (U, F) if £/ G L P (U,F) for any real-valued £ € C™{U). Also 
L p (£l) = R) and L p = L p ( R d ). We use to denote a definition. By T 

and F~ l we denote the d-dimensional Fourier transform and the inverse Fourier 
transform respectively, i.e. 


:= \ d /2 [ e lX< fi. x )dx , T 1 (f){ x ) 

(2ir) d / 2 J R d 


(27r) d / 2 




R d 


md£. 


For a Lebesgue set A C R d , we use |A| to denote its Lebesgue measure and by 
1,4(2;) we denote the indicator of A. For a complex number z, 9?[z] is the real part 
of z. Finally if we write N = N(a,b ,...), this means that the constant N depends 
only on a, b,.... 


2 . Main results 


We fix T G (0,00) throughout the article. For a > 0 denote 
ka{t) :=f x - 1 T{a)- 1 , t> 0, 

where T(a) := / 0 °° t a ~ 1 e~ t dt. For functions ip G Li((0, T)) the Riemann-Liouvillc 
fractional integral of the order a > 0 is defined as 

I a ip(t) = k a * <p(t) = —5— f (t- s) a ~ 1 (p(s)ds. 

r («) J 0 

It is easy to check that 

= I a l\{t) Va, /3 > 0. 


Also, by Jensen’s inequality, for any p G [ 1 , 00], 

ll-^VIlipdo.T)) < JV(r } a)|^|U p ((o ir )). (2.1) 

It is also known (see e.g. [39] ) that I a : B x —> C x+a is a bounded operator if A > 0 
and A + a < 1, where B° = Loo and B x = C x if A > 0. 





4 


ILDOO KIM, KYEONG-HUN KIM, AND SUNGBIN LIM 


Let fceN, k— 1 < a < k, and /( fe-1 )(t) be absolutely continuous, where /( fc_1 7) 
denotes the (k— l)-th derivative of function /. Then the Caputo fractional derivative 
of order a > 0 is defined as 


$7(t) = v(y ' /"*(* - s) fc -“-7 {n \s)ds 
T(k - a) J Q 


1 


(t - |/ (fc- 1 ) (s) - / (fe_ 1 ) ( 0 ) 


ds. 


( 2 . 2 ) 

(2.3) 


T(k — a) dt J Q 

Note that (12.311 (“Kochubei extension”) is defined for a broader class of functions. 
Let q > 1. For functions / £ C fc ([0,T]), we denote 

\ l / q / \ 1 / g 




l/W 








The following lemma shows that it is irreverent whether one uses (12.21) or m as 
the definition of <9“ for functions in H°(T). 


Lemma 2.1. The closures Hf(T) and H q (T) of C k ([0,T]) in the space L g ((0,T)) 
with respect to norms || • ||#c« related to (12.21) and (12.31) respectively coincide. 


Proof. This is obvious because ()2.2I) and (12.3[) are equal for functions / £ C fc ([0, T]). 

□ 


Next we introduce an another fractional derivative. Let Df denote the Riemann- 
Liouville fractional derivative of order a which is defined as 

D?<P(t) = |(J 1 -»(t), a £(0,1). (2.4) 

It is obvious that 

d?<p(t) = D?(ip - y(0)) = , a £(0,1). (2.5) 

It is easy to check for any ip £ Li((0,T)), 

D?I a <p = <p, a £(0,1). (2.6) 

Similarly, the equality 

I a D?<p = v, a £(0,1) (2.7) 

also holds if is absolutely continuous and I l ~ a ip{ 0 ) = 0 . 

Definition 2.2. Let k — 1 < a < k and f £ H q (T). We write /(0) = 0 if there 
exists a sequence f n £ C k ([0,T]) such that f n ( 0) = 0 and f n —> f in Hf(T). 
Similarly if a > 1 we write f'( 0) = 0 if f n { 0) = 0 for all n. 

The following lemma gives sufficient and necessary conditions for f £ Hf(T) 
and /( 0 ) = 0 (or /'( 0 ) = 0 ). 

Lemma 2.3. (i) Let a £ (0,1) and q > 1. Then f £ H q (T) and /(0) = 0 if and 
only if f £ L g ((0,T)), / 1_ “/ £ H q (T), I 1 ~ a f(t) is continuous, and I 1 ~ a f( 0) = 0. 

(ii) Let a £ (1,2) and q > 1. ITien / £ Hf(T ) and /'(0) = 0 if and only if 
f £ Hq(T), I 2 ~ a f £ H q (T), I 2 ~ a f'(t) is continuous, and I 2 ~ a f( 0) = 0. 
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Proof, (i) Suppose / £ H q (T ) and /(0) = 0. Take a sequence f n £ C' 1 ([0,T]) sat¬ 
isfying /„(0) = 0 and /” ->■ / in H*{T). Then we have /„->■/ and ^(/ 1 _Q 7„) -► 
9“/ in L g ((0, T)). It follows that I 1 ~ a f n converges to I 1 ~ a f in the space R 2 ((0, T)). 
It also follows that 

I l ~ a f(t) = f dtf(s)ds, t < T (a.e.). 

Jo 

Since 1 — 1/g > 0, by the Sobolev embedding theorem I 1 ~ a f(t) is continuous in t 
and thus the above equality holds for all t, which implies I 1 _Q /( 0 ) = 0 . 

Next we prove the “if’ part of (i). By the assumption, we can choose a function 
g £ L q ((0,T)) so that 

i 1 ~ a f(t)= f g(s)ds, vt<T. 

Jo 

Take a sequence of functions g n £ C 1 ([0,T]) which converges to g in L q ((0,T)). 
Define f n = I a g n . Then f n £ C^O,! 1 ]), /„(0) = 0, and 

fn^I a g = f in i,(( 0 ,r)). 

Also 

d?f n = d?I a g n =g n ^d?f in L q ((0,T)). 

(ii) The proof is very similar to (i). We only explain how one can choose a 
sequence to prove the “if’ part. 

By the assumption / £ H q (T) and q > 1, we may assume / £ C([0, T]). Take g £ 
L q ((0,T)) so that I 2 ~ a f = f*g(s)ds. We choose g n £ C 2 ([0, T]) which converges 
to g in L q ((0, T)). Define f n (t ) = I a g n (t ) + /(0). Then /£(0) = I a ~ 1 g n ( 0) = 0, 
fn —> f in L q ((0,T)), and dff n is a Cauchy sequence in L g ((0,T)). Thus / £ 
Hp(T) and /'(0) = 0. 

The lemma is proved. □ 


Next we introduce our solution space M q ’ k (T) and related notation. Roughly 
speaking, we write u £ H“’p (T) if and only if 

u, d t a u, D k x u £ L q ((0,T),L p ). 

For p, q > 1 and fc = 0,1, 2, • • •, we denote 

H k p = H k (K d ) ={ue L p (R d ) : Dfu £ L p (R d ), | 7 | < k}, 

Kp( T ) = L q ((0,T),HZ ), L,, P (T) := H 

where Df are derivatives in the distributional sense. Thus u £ H°’p(T) if and only 
if u(t,-) is H k -valued measurable function satisfying 


I u IIh“;£(t) := 


-I 1/9 


Ml H*ds 


< oo. 


We extend the real-valued time fractional Sobolev space to i p (R d )-valued one. In 
other words, we consider the completion of C 2 ([0,T] x R d ) nL w (T) with respect 
to norm 

II ' IIl 9iP (T) + II &t ' IIl, iP (T) 

in Lq, p (T). 
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Definition 2.4. For a £ (0,2) we say u £ H“’p(T) if and only if there exists a 
sequence u n £ C 2 ([0,T] x R d ) flL w (T) so that sup n ||<9f tt n |k,, (t) < oo, an d 

\W - Wnlk,,„(T) -t 0 and ||9fu n “ 5 t Q U m || L ,, p (T) 0 

os n and m qo to infinity. We call this sequence u„ a defininq sequence of u. For 
u £ H“;°(T), we define 

<9“it = lim <9“u n in L op (T), 

n—>• oo 

where u n is a defining sequence of u. Obviously 11“’°(T) is a Banach space with 
the norm 

II u IIh“; p °(t) = IMk„, P (T) + II^“w||l,, p (t)- 

Definition 2.5. for u £ H“’°(T), we say that it(0, x) = 0 if any only if there exists 
a defining sequence u n such that 

u n (0, x) = 0 Vx £ R d , Vn £ N. 

Similarly we say that u(Q, •) = 0 and Jjit(0,-) = 0 if any only if there exists a 
defining sequence u n such that 

d 

u n ( 0,x)=0 and —u n (0,x) = 0 Vx £ R d , Vn £ N. 

Let 

H“; p fc (T) := H“; p °(T) n E° q k p (T) 

and H“’p |0 (T) be the subspace of (T) such that 

u(0, •) = 0 if a £ (0,1] 

and 

it(0,-) = 0 and J^u(0, •) = 0 if a £(1,2). 


Theorem 2.6. (i) The space H(T) is a Banach space with the norm 

IMIiH“’ p fc (T) := H u IIh°; p (t) + II m IIh“’°(t)- 


(ii) The space H“p 0 (T) is a closed subspace o/H“’p(T). 

(Hi) C'£°(R(J_ +1 ) is dense inH“p 0 (T). 

(iv) For any u £ H“ p 2 0 (T) ; 

\\u(t)\\ Lp <N(a) [ (t- s)“" 1 ||5 i “n(s)|| Lp ds, t<T (a.e.). (2.8) 

Jo 

Consequently, 


\\u\\l qAt) <N{q,a,T) j (s - r) a 1 \\d?u{r)\\ q Lp drds, \/t<T. 

Proof, (i) This is obvious because both H“’°(T) and H°^;(T) are Banach spaces. 

(ii) Suppose u n £ H“'p 0 (T) and u £ H“’p (T) so that u n — > u in H“’p (T). Since 
u n £ 0 (T) for each n, we can find a v n £ C 2 ([0, T] x R d ) n L 9 jP (T) so that 

IK - Unlk ? , P (T) < - and ||9“n„ - dKlk,, p( T) < K 


v„(0,x) = 0 Vx £ R rf , if a £ (0,1], 
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and 

d 

u„(0,x) = 0 and —v n (0,x)=0 Vx € R d , if a €(1,2). 
Therefore u £ H“’^ \ 0 (T) because obviously 

|| V n — U 11 jja p (T) ~t 0 aS 12 00 ■ 

This certainly proves (ii). 

(iii) We take nonnegative smooth functions 77i £ C'^°((1, 2)), 772 £ C£°(R d ), and 
773 £ C~(R d ) so that 

/ 771 (i) dt = 1 , / 772(re) dx = 1 , and 773(2;) = 1 if |x| < 1 . 

20 Jr j 

For £ 1 , 62,83 > 0, we define 


Vl,eAt) =£l 1 m(i/ £ l)» m,e(x) = £ 2 d ri 2 {x/£ 2 ), 

/•OO 

u El (t,x) = / u(s, x)rji tEl (t — s)ds, 

Jo 

UE2 (t,x) = / u(s,y)r) 1>ei (t - s)r] 2 , £2 {x - y)dsdy, 

J R d 20 


and 


u El ’ E2 ’ Es (t, x) = 77(^)773 (£ 32 ) / / u(s, y)rj 1<Sl (t — s)r] 2 ,e 2 (x — y)dsdy, 

Jn d Jo 

where 77 £ C'°°([0,oo)) such that 77 (f) = 1 for allt < T and vanishes for all large t. 
Due to the condition 771 £ C£°((l, 2)), it holds that 

u £l ’ e2 ’ e3 (t,x) = 0 Vt < £1, Vx £ R d . 

We can easily check that for any u £ H“’^ 0 (T) 

d?u E At) = (d?u) ei (t). 

Hence for any given e > 0 we have 


||u — u El ’ E2 ’ £3 |l H “’ p fe (T) 

— II U — uEl ll]H“’ p fe (T) + H^ 61 — uei ’ e2 HlH“’p (T) + ||w ei ’ £2 ^ U ei,E2 ’ E3 \\ m *,k ( T ) < £ 

if £1, £2, and £3 are small enough. Therefore (iii) is proved. 

(iv) Due to (iii), it is enough to prove (12.811 only for u £ C)? 0 (R d+1 ). Denote 
/ := dfu. One can easily check 


u(t) = f k a {t — s)f(s)ds, \/t<T 
Jo 

in the space L p , which clearly implies (12.811 due to the generalized Minkowski in¬ 
equality. The theorem is proved. □ 


Assumption 2.7. Let f(u ) = f(t,x,u) and f 0 = f(t,x,0). 

(i) There exist 0 = To < T) < ■ ■ ■ < Tf = T and functions ( t, x) such that 

i 

a lJ (t,x) = J 2 a%(t,x) I {Tk _ u T k ](t), (a.e.). 
k =1 
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(ii) There exist constants 6 , K > 0 so that for any k, t, and x 

m 2 <m 2 , V£eR d , (2.9) 

and 

\ a k(t,x)\ + \b l ft,x)\ + \c(t,x)\ < K. 

(Hi) The coefficients a^ are uniformly continuous on ftk-i,tk) x R d for all k = 
1 ,...,t and i,j = 1 ,... ,d. 

(iv) /o £ L q, P {T) and ffa) satisfies the following continuity property: for any e > 0, 
there exists a constant K e > 0 such that 

II f(t,x,u) - f(t,x,v)\\ Lp < e\\u-v \\ H 2 + K s \\u-v\\ Lp , (2.10) 

for any ft, x) and u, v £ H 2 . 

If p ^ q then we need an additional condition (see the comment below (15.411 for 
the reason). 

Assumption 2.8. If p ^ q then lim| ;i; |_ i . 00 a lJ (t,:r) exists uniformly in t £ (0, T). 

Here is the main result of this article. The proof will be given in Section [5] 

Theorem 2.9. Letp,q > 1 and Assumptions [377] and HTSI hold. Then the equation 

dfu = a^u x i x j + b l u x i + cu + f(u), t > 0 (2-11) 

admits a unique solution u in the class H q'p 0 {T), and for this solution it holds that 

IMIh“’p(t) — ^o||/o||l ? , p (t)) (2.12) 

where Nq depends only on d,p , q , <5, K , K e ,T , £, and the modulus of continuity of . 

Remark 2.10. (i) Due to the definition of our solution space H“p 0 (T), the zero 
initial condition is given to equation (12.111) . that is 

du 

u( 0 ,x) = 0 and additionally — (0,x) =0 if a > 1. 

(ii) Some examples of f(u) satisfying (12.101) can be found e.g. in [19]. For 
instance, let n := 2 — d/p > 0, h = h{x) £ L p: and 

f(x,u) := h{x) sup |u|. 

X 

Then by the Sobolev embedding theorem, 

\\f(u) - f(v)\\ Lp < \\h\\L p sup \u - v\ < N\\u-v\\ h * < e\\u - v \\ H 2 + K\\u - v\\ Lp . 

X 

Similarly one can show that f(t, x, u) := aft, x)(—A) s u also satisfies (12.101) if aft, x) 
is bounded and <5 £ (0,1). 
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3 . Some Preliminaries 


In this section we introduce some estimates of kernels related to the representa¬ 
tion of a solution, and we also present an /^-estimate of a solution. 

The Mittag-Leffler function E a (z) is defined as 


E a {z) = ]T 


fc=o 


T(afc + 1) 


, z £ C, a > 0. 


The series converges for any z £ C, and E a (z) is an entire function. Using 

,/3 _ r(/3 + 1) 


f¥*fP = _ ~ _ +P~ 

‘ T(/3 + 1-a) 


ft > a 


one can easily check that for any constant A, 


(pit) := E a (\t a ) 

satisfies <^(0) = 1 (also y/(0) = 0 if a > 1) and 


dfip = A ip, t > 0. 


By taking the Fourier transform to the equation 

dfu = Ait, t > 0, u(0) = h, (andit'(O) = 0 if a > 1) 

one can formally get E(u)(t,f) = E a (—t a \f\ 2 )h. Thus it is naturally needed to find 
an integrable function p(t,x) satisfying JF(p{t, •))(£) = E a (—t a \f\ 2 ). It is known 
that (see e.g. (12) in [TUI Theorem 1.3-4]) 


Ea( ~ t] ^ iT\tV t>0 ' (3 - 1} 

Therefore, E a (—t a |£| 2 ) is integrable only if d = 1. If d > 2 then tF~ l {E a {— i“|£| 2 )) 
might be understood as an improper integral. However, in this article we do not 
consider F~ 1 (E a {—t a |£| 2 )). 

Lemma 3.1. (i) Let d > 1 and a £ (0,2). Then there exists a function p(t,x) 
such that p(t , •) is integrable in R d and 

T(p(t,-m = E a (-t a ici 2 ). 

(ii) Let m,n = 0, 1,2, and denote R = f -a |:r| 2 . Then there exist constants C 
and a depending only on m, n, d, and a so that if R > 1 

\d?D?p(t,x)\ < Af^^-"exp{_ar^|*|^}, (3.2) 

and if R < 1 

| d?D2p(t,x)\ < N\x\- d - m t~ n (R + R\nR- l d=2 , m=0 + R 1/2 • l d =i, ro=0 ) . (3.3) 


By (13.21) . p{t , x) is absolutely continuous on (0, T) and lim^o pit, x) = 0 if x ^ 0. 
Thus we can define 


q(t,x) 


I a ~ l p{Lx ), a £( 1 , 2 ) 
D\~ a p{t,x), a £(0,1). 


Since p(0, x) = 0 for x 0, D\ a p(t, x) = d] a p(t , x). 
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Lemma 3.2. (i) Let d > 1, a £ (0, 2), and m, n = 0,1, 2, • • •. Denote R = t a |x| 2 . 
Then there exist constants N and a depending only on m, n, d, and a so that if 
R> 1 

\d?D™q(t,x)\ < Nt ~“ ( " +m) -"+<»- 1 exp{-trt~^ \x\ ^}, (3.4) 

and if R < 1 

| d?D™q(t, a:)| < N\x\- d - m t- n+a -\R 2 + R 2 Ini? • l d=2 ) 

+ N\x\~ d t~ n+a ~ 1 (r 1/2 ■ l d= i + R ■ l rf=2 + R 2 In R ■ l d=4 ) 1 m=0 . 

(3.5) 

(ii) For any t ^ 0 and x ^ 0, 

= Ap, ^ = Ag. (3.6) 


One can find similar statements of Lemmas 13.11 and 13.21 in [TTJ [121 HU 133 ■ For 
the sake of completeness, we give an independent and rigorous proof in Section [5] 


Corollary 3.3. Let 0 < e < T. Then 


and 


/ sup \D™q(t, x)\dx < oo, m = 0,1,2. 
JR d te[e,T] 


/ sup 

m d te[e,r] 


| D™p(t, x)\dx < oo, 


m = 0,1. 


Proof. By (13.41) and (13.51) . it follows that for large |x| 

sup \D?q(t,x)\ < Ne~ c W 

t€.[e,T] 


and for small |x| 

sup \D™q(t,x)\ < N\x\~ d ~ m (|x| 3 + |x| • l m=0 ), 

t£[e,T] 

where N and c depend only on d, a, T, and e. This certainly proves the assertion 
related to q , and p is handled similarly. The corollary is proved. □ 


Lemma 3.4. Let A > 0, a € (0,2), and if be a continuous function on [0,oo) so 
that the Laplace transforms of if and d’fcf exist, 

d?ct> + \cj> = f{t), t> 0, (3.7) 

and <(>(0) = 0 ( additionally <//(0) = 0 if a £ (1,2)). Moreover we assume for each 
s > 0, 

e~ st f df<p(r)dr —> 0 as t —> oo. 

Jo 

Then 

<f(t) = [ s)f(s)ds, (3.8) 

Jo 

where H a ,\(t) = D]~ a E a (—Xt a ) if a £ (0,1) and H a ,\{t) = I a ~ 1 E a (—Xt a ) other¬ 


wise. 
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Proof. First, recall that ip(t) := E a {—t a ) satisfies </?(0) = 1 and dfip = —<p. Let C 
denote the Laplace transform. Then from C[d^^p] = —£[</?] we get 

roo „a-l 

Jo s + i 

Here, we used the following facts: for fi £ (0,1), 

C[h'}= sC[h](s), C[h* g](s) = C[h\(s) ■ C[g}(s), £[fci-/s](s) = s' 3-1 . (3.9) 

It follows that 

C[E a (-\t a )](s) = C[H a , x ](s) = -J—. (3.10) 

S t A S ~r A 

On the other hand, taking the Laplace transform to (13.71) and using (13.91) we get 

£[<£](«) = * x ' £[/](«), s > o. 

This and (13.101) certainly prove the lemma, because to prove equality (] 3.8 D it is 
enough to show that two functions under consideration have the same Laplace 
transform. □ 


Lemma 3.5. (i) Let u £ C£°(R+ +1 ) and denote f := d?u — A u. Then 

u(t, x)= / q(t-s,x-y)f(s,y)dyds. (3.11) 

Jo J R d 

(ii) Let f £ C£°(R+ +1 ) and define u as in (13.111) . Then u satisfies <9“m = A u + f. 
Proof, (i) Since q is integrable on (0, T) x R d (see Lemma EPl) . it is enough to prove 

u(t,£)=[ q(t- s,€)f(s,£)ds, (3.12) 

Jo 

where / denotes the Fourier transform of / with respect to x, i.e. 

/(s, £) := E(f(s, ■))(£) = , \ d/2 f e~ lx< f(s, x)dx. 

First note that from the definition of / we get 

f{t, 0 = d?u(t,£) + |f| 2 u(t,C), > °- 

Therefore by Lemma 13.41 

u{t,Q=[ H at \tf(t-s)f{s,€)ds. 

Jo 

Hence it is enough to prove 

q(t,0=H a ^2(t). (3.13) 

Denote Cd = (27r) _d / 2 . If a £ (0,1) then by the definition 


q(t,x) = Cd [ e lx '^q(t 1 x)dx = Cd [ e J- [ k a (t - s)p{s,x)ds 
Jn d JR d i at J o 


dx. 
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Since q(t, •) is integrable in R d uniformly in a neighborhood of t > 0 (see Corollary 
13.31) . one can take the derivative ^ out of the integral. After this using Fubini’s 
theorem we get 


Q = 


dt 


k a (t. s) 


Cd 




/R d 


p(s, x)dx 


ds = D\- a E a (—t a \£\ 2 ). 


Hence (13.131) and (i) are proved. The case a £ [1, 2) is easier and we skip the proof. 


(ii) Taking the Fourier transform to (13.111) we get 

u{t,i)= ( H a ^(t-s)f(s,g)ds. 

Jo 

Note that if one defines 0 as in (|3.8[) then it satisfies (|3.7f) . Consequently, u satisfies 

+ |£| 2 u{t,C) = /(*,£)> 

and this certainly proves the equality d^u = A u + f, because f(t, •) = 0 if t is small 
enough and thus for each t > 0, 

F (d?u(t, •)) (x) = [ H a> \z\ 2 {t- 8 )d?f{s,t)ds = d?u(t,£). 

Jo 

The lemma is proved. □ 


Now we define the operator Q by 


Gf{t,x) 



s,x- y)Af(s,y)dyds. 


Since q is integrable on (0,T) x R d , Qf is well-defined if / £ C(i 0 (R rf+1 ). Also 
recall that D x q{t : ■) and D x q(t, •) are integrable in R d for each t > 0, and therefore 
it follows that 


Gf(t,x) = lim 

e —>0 




Aq(t -s,x- y)f{s,y)dy 


ds. 


(3.14) 


Lemma 3.6. Let f £ C)? 0 (R' i+1 ). Then 

II£/IIl 2( r^) < tf||/IU a( R«+i) (3-15) 

where N = N(a,d). 

Proof. Denote qM = <?lo<t<M- One can easily check that qM is integrable in R d+1 . 
Denote Gm f = qM * A/. Then by Parseval’s identity 

\\G M f\\l 2 = f \T d+1 (q M *Af)\ 2 drdC 

dR<G-l 

By the properties of the Fourier transform, 

T d+1 (q M *Af)(T,Z) = -N(dm 2 T d+1 (q M )(p^d + i(f)(r,0, V(r,£) £ R d+1 . 

Set 

:= -|^| 2 ^d+i(?M)(r,0 = [ e~ lTt \F d (^)(t^)dt. 

J 0<t<M ot 

Now we claim that /m(dO is bounded uniformly for M > 0, i.e. 

sup |/m(t,£)I < 00- 

M>0,r,£ 
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Then the claim implies 

\\GMf\\l 2 <N\\f\\l 2 , (3.16) 

where N is independent of M. 

By the integration by parts and the change of variables, 

r-rM i 

Jm(t, 0 = *sgn(r) / 0 )“||£| 2 )dt + e~ iTM E a (-M a \^\ 2 ) - 1. 

Jo \ T \ 

(3.17) 


Denote for ^ < 77 < min{ 7 r, cnr} 

A* := {z £ C : | 7 r — Arg 2 r| <71 — 77 }. 

Then by [TO] (1.3.12)] or 01) Theorem 1.6], 

|E ” WI = TTN’ zeA; ' 

Hence the function E a (z) is bounded in A*. 

Let r > 0. We take a 9q £ (0, 7t/ 3) sufficiently small so that —e~ q,e , —e ~ la9 £ A* 
for any 9 £ [0,$o]- Denote 

C'i.rM = {t : t £ [0, tM}}, C 2 ,tM = {te~ l9 ° : t £ [0, tM}}, 


C 3 , t m = {rMe~ i0 : 9 £ [O,0 o ]}, 

and define a contour C t m = Ci,rM U C 2 , t m U C 3 , t m- Then the contour integral 
of the function e~ lz E a (— (y)“|£| 2 ) on C t m is zero for any M > 0. Since E a (z) is 
bounded in A*, 


>c 3 ,- 


e-'*E a (-{-) a \tf)dz 

T 


< 


(rM)e 


—tM sin 6 


de < n( 0 o ). 


(3.18) 


Also, 


e- iz EJ-(-r\tf)dz 


prM 

<N g-tsin e 0 dt < n{6 0 ). 


(3.19) 


JC2,tM t do 

Note that 9 0 depends only on a. Hence by (13.171) . (13.181) . and (13.191) . it follows that 
if r > 0 then \Im(t, £)] is bounded uniformly for M. 

If r < 0 then we choose a contour C’ tM = Ci jT m U C 2 tM U C' 3 tM where 

C 2 .rM = {te l9 ° : t £ [0, tM]}, C' >tM = {rMe ie : 6 £ [0,0 O ]}- 

Then the same arguments above go through. Thus our claim is proved. 

To finish the proof, observe that for each (f, x) £ R d+1 


q M * A/(t,x) = g M f{t,x) —>• Gf(t,x) as M-» 00 
since A/ £ C^°(R d+1 ). Therefore (13.161) and Fatou’s lemma easily yields (|3.15l) . □ 

















14 


ILDOO KIM, KYEONG-HUN KIM, AND SUNGBIN LIM 


4. BMO AND L g (L g )-ESTIMATE 

In this section we obtain BMO and L g (L p )-estimates for a solution of the model 
equation 

dfU = A u + /, (t, x) £ (0, oo) x R d . 

Recall that p is integrable with respect to x, E{p{t, •)) = E a (—t a |£| 2 ), and q is 
defined as 

I a ~ l p{t,x ), a £ ( 1 , 2 ) 

D\~ a p{t,x), 0 6 ( 0 , 1 ). 


ds. 


(4.1) 


q(t,x) = 

Also recall that for / £ C)? 0 (R t | +1 ), it holds that 

Gf = [ [ A<?(f- s,x-y)f(s,y)dy 

J — oo 

and by (|B. 15|) the operator Q is continuously extended onto L 2 (R d+1 ). We denote 
this extension by the same notation Q. 

For a locally integrable function h on R d+1 , we define the BMO semi-norm of h 
on R d+1 as 

II^IIbmo(R‘ 1 + 1 ) = SU P T 7 m / \ h (t, x ) - h Q \dtdx 
QeQ \Q\ JQ 

where fiQ = h(t,x)dtdx and 

Q := {Qs(t 0 ,x 0 ) = (to ~ S 2/a ,t 0 + 6 2/a ) x B s (x 0 ) : S > 0 , (t 0 , x 0 ) £ R d+1 }. 
Denote Qs ■= Qs( 0 , 0 ). 

Lemma 4.1. Let f £ L 2 (R d+1 ) and vanish on R d+1 \ Q 3 S- Then 

/ \Gf(t,x)\dtdx < 

JQs 

where N = N (d, a). 

Proof. By Holder’s inequality and Lemma 13.61 

[ \Gf(t,x)\dtdx < ( f \Gf{t, x)\ 2 dtdx\ |Qa| 1/2 

JQs \Jq 6 / 

< \\Gf\\ L 2 i n^)\Qs \ 1/2 < N\\f\\ L 2 (n^)\Qs \ 1/2 
= N [J q 1/(^1 x)\ 2 dtdx^j \Qs\ 1/2 < N\\f\\ LaoiRd+1) \Qs\. 


The lemma is proved. 


□ 


Denote 

K(t, x) = lt > 0 Ag(t, x). 

Due to Lemma [3A)ii), we have K(t,x) = lt>o d t p(t,x). Furthermore the following 
scaling properties hold (see (16.2111 and (16.121) for detail): 

K(t, x) = AT(1, t~ a / 2 x) (4.2) 


d t K(t,x) = l t>0 t 2 a * d (d t K)(l,t a/2 x), 


( 4 . 3 ) 
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and 

d 1 c*(d+l) d /o 

g-K{t,x) = l t >ot 2 g-K(l,t a / x). 

Lemma 4.2. There exists a constant N = N(a,d) such that 
(i) for any t > a and rj > 0, 



s, y)\dyds < N(t — a)rj 2 /“; 


(ii) for any t > r > a, 


r / 

J-o o JR d 


- s,y)~ K(t 


s, y)|dyds < iV 



(4.4) 


(4.5) 


(4.6) 


- 2 


(m) /or any t > a and x 6 R d , 

f f \K(t — s,x + y) — K(t — s,y)\dyds < N\x\(t — a) 

J -oo JR d 

Proof. First observe that 

(|y| 2 / Q |7F(l,y)| + |^(1,y)| + |d t /v (1, y)|) dy < oo 

1-1 \ ' 


(4.7) 


/Rd+i 


which is an easy consequence of Lemma 13.21 By (14.21) , it holds that 

f f \K(t-s,y)\dyds= f (i-s) -1- ^ f K(l, (t - s)~ a/2 y)dyds 
Ja J\y\>n Ja J\y\>n 


\y\>v 

= f (t-s)- 1 [ K(l,y)dyds 

Ja J\y\> r i(t—s )~ a / 2 

r t 


< 


J V~ 2/a ds ) (J \y\ 2/a K{l,y)dy 


< N(t-a)ir 2/a - 

Hence (i) is proved. 

Next we prove (ii) and (iii) on the basis of the scaling property. By (14.31) , (14.41) , 
and the mean-value theorem, we have 

[ ( I K{t-s,y)~ K{t - s,y)\dyds 

J- oo J R d 

< (t — t) f f f \d s K(8t + (1 — 9)t — s, y)\dddyds 

J —oo J R d J 0 

(/ ( dt + ( 1 ^ d ) T ^ s )~ 2ds ^j (^J d \ d sK{l,y)\diJj dO 

(/ (r — s)~ 2 ds S j < N-^ 




t — r 
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Also, 


[ [ \K(t — s,y + x) — K(t — s, y)\dyds 

J — oo J 


< \x\ f f f \D y K(t — s,y + Qx)\d6dyds 
J — oo J R d J 0 

<\x\ f f \D y K(s, y)\dyds 

Jt-a JR d 

<\x\ [ [ |lt>oi -1- ^-DyAT(l, y))\dyds < N\x\(t — a) 

J t—a J 

The lemma is proved. 


_ «W “/ 2 


Lemma 4.3. Let f £ L 2 (R d+1 ) and f = 0 on Q 25 - Then 

/ / \Gf(t,x) - Gf(s,y)\dsdydtdx < N{d,a)\\f\\ Laa(:R d+ 1 y 
JQs'' Qs 

Proof. First assume / £ C'“(R d+1 ). We claim that 

/ | Sf(t,x) - Gf(-S 2/a ,0)\dtdx < 7V||/|| ioo(R , +1) . 

JQs 

Let (f, x) £ Qs- Then 

< | Sf(t, x) - Qf(t, 0)1 + I gf(t, 0) - Gf(-S 2/a , 0)1 =: 1 1 + J 2 . 
We consider first. 

rt 


h = 


[ ( s,x-y)~ K(t-s,-y))f(s,y)dyds 

J- oo J RA 


r r-( 25) 2/o1 i- 

/ •• • dyds+ / • • • dyds 

/-(2<5) 2 /“ J R d J -oo JR d 

< [ [ \K(t-s,y)\\f{s,x-y)\dyds 

2<5) 2 /“ J R d 

+ [ [ \K(t- S,y)\\f(s, -y)\dyds 

7 _( 25 ) 2 /c J R d 
2 S) 2/a r 

+ / \K{t- s,x- y) - K(t- s, -y)\ \f(s,y)\dyds 

J —oo J 

=: In +112 + ^13- 

Note that if —(2<5) 2 /“ < s < t < 8 2 / a and \y\ < 5, then 

f{s,x-y)= 0, f(s,-y) = 0, 


□ 


(4.8) 


(4.9) 


(4.10) 
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because \x — y\ < 25, | — y\ < 5, and / = 0 on Q 25 ■ Hence by (14.51) . In +I 12 is less 
than or equal to 

N \\f\\L oc ,(R d + d ) f f \K(t — s,y)\dyds 

J-(2S) 2 /<* J\y\>6 

< JV||/|| iao(R< i + i)(t + (2<5) 2 /“)<T 2 /“ < iV||/|| ioo(Rd+1) . 

Also, by (ECU) we have 

/• — (2<5) 2/ “ r 

X 13 < ^V||/||i oo(R d+i ) / / \K(t-s,x-y)~ K(t-s,-y)\dyds 

J — 00 J R d 

< JV||/|| iao (R-+i)|*|(t + (2<5) 2 /“)-“/ 2 < 7V||/|| ioo(R , +1) . 

Next, we consider I 2 . Note that 


I 2 = 


< 


+ 


Gf(t,0)-gf(-5 2 ' a ,0) 

[ [ K{t-s,y)f(s,-y)dyds- [ [ K(-5 2/a - s,-y)f(s,-y)dyds 

J- 00 J R d J- 00 J R d 

t p p —<5 2 / a p 

/ K(t-s,y)f(s,-y)dsdy- / K(t - s,y)f(s,-y)dyds 

JR d J- 00 JR d 




— OO 

—< 5 2 / 


/R d L 


K(t — s, y) - K(-5 2/a - s, y)j f(s,-y)dyds 


=: X21 + 122- 


Recall (14.101) . Thus by (14.51) . we have 


I21 < 


-5 2 /“ -/R d 


|AT(i- s,y)/(s,-y)|dyds 


Also, 


I22 < 


< ll/llL 00 (R d + 1 ) [ [ \K(t-s,y)\dyds 

J-S 2/“ I|y|>5 

< AT||/|| ioo(Rd+1) (t + <5 2/ “)<5- 2 /“ < A'||/|| Loo(R d+ 1) . 

_ ^2/a 

[ [ \ K (t~ s,y)~ K{-5 2/a - s,y)\\f{s,-y)\dyds 

J-(26) 2 / a JR d 

r-( 2 S) 2/a r 

+ ll/llWR-+ 1 ) / / K{t~s,y)~K{-5 2 / a -s,y) 

J— 00 «/R d 


dyds 


—: X221 +1222 • 
By (14.51) again, we have 

Z221 < ll/ll. 


(R d +!) / / |AT(t - s,t)|dyds 

J-( 2 < 5 ) 2 /“ J|t |><5 

+ ll/llL 00 (R d +i) [ [ \K{-S 2/a - s,y)\dyds 

J-(25) 2 /°‘ J\v\>6 


< N\\f\\ LUR ^y 
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On the other hand, by (14.6[) we obtain 

t + (5 2 / a 

1222 - _J2/a + (2^)2/a ll/ll^oo(R d + 1 ) < ^ll/lll^R^ 1 )' 

Hence (14.91) is proved and this obviously implies (14.81) for / £ C'^ 0 (R d+1 ). 

Now we consider the general case, that is / £ L 2 (R d+1 ). We choose a sequence 
of functions f n £ C'“(R d+1 ) such that /„ = 0 on Q 28 , Gfn —t Gf (a.e.), and 
ll/nllwRrf+i) < II/IILoo(R d + X ) • Then by Fatou’s lemma, 

/ / |£/(*,a:) - f//(s,y)| dtdsdxdy 
JQs J Qs 

< liminf f f \Gf n (t, x) - Gf n (s, y)\dtdsdxdy 
n ^°° JqJqs 

< IVliminf \\f n \\L ao (R d+1 ) ^ JV H/IU= 0 (R d + 1 )- 


The lemma is proved. 


□ 


(4.14) 


Theorem 4.4. (i) For any f £ L 2 (R d+1 ) D I/ 00 (R d+1 ), 

ll^/llsMO(R d + 1 ) < IV (d, Cl) ||/||/, 00 (R d + 1 ) • 

(ii) For anyp,q£ (l,oo) and f £ C)? 0 (R d+1 ), 

II^/IIl,(r.l p ) < N (d,p,q,a)\\f\\ Lq ( RLp y (4-12) 

Proof, (i) It suffices to prove that for each Q = Qs{to,Xo) the following holds : 

I Gf{t,x) - {Gf) Q \dtdx < A r ||/|U 00 (R<i+i). (4.13) 


/ 

Jq 


Due to the translation invariant property of the operator Q , we may assume that 
(to, Xq) = 0. Thus 

Q = Q s = (-S 2 / a ,S 2/a )x B S . 

Take £ £ C^°(R d+1 ) such that ( = 1 on Q 25 and ( = 0 outside Q 3 g. Then 


/ I Gf(t,x) - [Gf) Q \dtdx 
Jq 

< / \G(Cf) ~ (<7(C/))q | dtdx + / |0((1 - Of) ~ m 1 - C)/))q I dtdx 
Jq Jq 


'Q ->Q 

< 2 / \G((f)\dtdx + -f l |<?((1 - Of)(t,x) - £((1 - Qf){s,y)\dsdydtdx. 

Jq JqJq 

Thus (14.131) comes from Lemma l4~Tl and Lemma T4.31 

(ii) Step 1. We prove (14.121) for the case q = p. First we assume q = p > 2. For 
a measurable function h(t,x) on R d+1 , we define the maximal function 


Mh(t, x) = sup 4- \f(r,z)\drdz, 
q&qJq 


and the sharp function 


h^(t, x) = sup -f | f{r,z) - f Q | drdz. 
q&qJq 
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Then by the Fefferman-Stein theorem and the Hardy-Littlewood maximal theorem 
ll^ tt ||i ! 3 (R d + 1 ) ~ ll- A ^ /l llL p (R‘ i + 1 ) ~ ll^lli p (R‘ i + 1 )- (4-14) 

Combining Lemma [3761 with (14.141) . we get for any / £ L 2 (R d+1 ) D Loo(R d+1 ), 
IK^/)“llL 2 (R^)<iV||/|U 2 (R, +1 ). 

Also, by gUD 

UGfyh^n^) < N\\f\\ Loo{Ra+1) . 

Note that the map / —»■ ( Qf )** is subadditive since Q is a linear operator. Hence by 
a version of the Marcinkiewicz interpolation theorem (see e.g. PH Lemma 3.4]), 
for any p £ [2, oo) there exists a constant N such that 

II (^/)^ llip(R d + 1 ) < ^ll/llip(R d + 1 )> 

for all / £ L 2 (R d+1 ) LI Loo(R d+1 ). Finally by the Fefferman-Stein theorem, we get 

ll^/llip(R ci + 1 ) < N(d, a,p)||/||ip(R<i+ i ). 

Therefore (14.121) is proved for q = p £ [2, oo). 

For p £ (1, 2), we use the duality argument. Let /, g £ C'^°(R d+1 ) and p' be the 
conjugate of p, i.e. p' = £ (2,oo). By the integration by parts, the change of 

variable, and Fubini’s theorem, 


/Rd+i 


g(t, x)Qf(t , x)dxdt 


/ &g(t,x) ( / l t > s q(t — s, x — y)f{s , y)dsdy ] dxdt 

J R d +! \9R d + 1 / 

/ f(s,-y)[ l a>t q{s-t,y-x)Ag(-t, -x)dxdt) dsdy 

J R d +! \JR' 1 + 1 / 

f{~s,-y)Gg(s,y)dsdy 


/r ^ 1 


(4.15) 


where g(t,x) = g{—t, —x). Then by Holder’s inequality (also recall 2 < p' < oo), 


/Rd + i 


g(t , x)Gf{t , x)dxdt 


< mfhJGgh. < N\\f\\ Lp \\g\\ Lpl . 


Since g £ C')? 0 (R ti+1 ) is arbitrary, (14.111) is proved for <7 = p £ (1,2) as well. 

Step 2. Now we prove (|4.12l) for general p, q £ (1, 00 ). For each (i, s) £ R 2 , we 
define the operator /C(t, s) as follows: 

1C{t,s)f(x):= f K{t - s,x - y)f{y)dy, feC™. 

J R d 

Let p £ (1, 00 ) and (f, s) £ R 2 . Then by the mean-value theorem, Lemma 13.21 and 

(TO) . 


mt,s)f\\ Lp = 


f K(t~ s,x-y)f(y)dy 
JR d 


< 


\ —11 


Lp [ \K(t — s,y)\dy < N(d, a)(t — s) 

J R d 

Hence the operator /C(f, s ) is uniquely extendible to L p for t ^ s. Denote 
Q := [to, t 0 + <5), Q* := [to — <5, to + 25). 


Lt> * 
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Note that for t Q* and s, r £ Q, we have 

|s — r\ < S, 11 - (t 0 + (5)| > S, 

and recall K(t — s, x — y) = 0 if t < s. Thus by Lemma f3. 2 1 and (14.31) . it holds that 
\\IC(t,s)-IC(t,r)\\ HLp) 


= sup 


[ {K(t-s,x-y)-K(t-r,x-y)}f(y)dy 

J R d 


< sup \\f\\ Lp [ \K(t-s,x) - K(t-r,x)\dx 

imu P =i 

f i T s -/ v r w M7 N(d,a)\s — r\ 

= L mt ~ s ’ x) - K(t - r - x)ldx s «-fa + j)A 

where L (L p ) denotes the bounded linear operator norm on L p . Therefore, 

ls - r| dt 


^ {t _ (to + s) y. 


< iV|s — r\ f 

J\t 


1 


|t—(t 0 +5)|>5 (t _ ^0 + ^)) 2 


dt < NS 


r 2 dt < N. 


Furthermore, by following the proof of Theorem 1.1 of [20], one can easily check 
that for almost every t outside of the support of / € C'^ XD (R, L p ), 


/ OO 

IC(t,s)f(s,x)ds 

-OO 


where Q denotes the unique extension on L p (R rf+1 ) which is verified in Step 1. 
Hence by the Banach space-valued version of the Calderon-Zygmund theorem [201 
Theorem 4.1], our assertion is proved for 1 < q < p. 

For the remaining case that 1 < p < q < oo, we use the duality argument. 
Define p' = p/(p — 1) and q' = q/(q— 1). Since 1 < q' < p', by (14.151) and Holder’s 
inequality, 


/ R d +! 


g(t,x)Q f(t,x)dxdt = 


/R 


' R d 


f(-s,-y)Gg(s,y)dy I ds 


< [ \\f(-s,-)\\L p \\G~g(s,-)\\L p ,ds 

< N(d, a,p)||/||i 9 (R,L p ) ||g||L 9 ,(R,L p /) 
for any f,g£ C™( R d+1 ). Since g is arbitrary, we have 

||!?/I|l,(r,l p ) < N(d,a,p)\\f\\ Lq ^ Lp )- 

The theorem is proved. 


□ 


5. Proof of Theorem 12.91 

First we consider the solvability of the model equation. 

Lemma 5.1. Theorem \2.9\ holds for the equation <9“u = A u + fo with No = 
N 0 (p,q,a,T). 
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Proof. First we prove the uniqueness. Suppose that u £ H“ p 0 (T) and dfu = A u. 
Due to Lemma [Ml (iii), there exists a u n £ C^°(R d+1 ) such that u n —> u in 
i“p 0 (T). Due to Lemma [3~5l fib 


U n (t,x) 


S,x- y)fn(s,y)dyds , 


J o J R d 

where f n := dfu n — A u n . Since f n —> 0 in L|’(T), we obtain the uniqueness from 
Lemma 12.61 (iv) and Theorem 14.41 

For the existence and (12.121) . first assume / £ C£°(R+ +1 ) and define u = Qf. 
Then all the claims follows from Lemma ldAT iib Theorem l4.41 and (12.81) . For general 
/ one can consider an approximation /„—►/, and above arguments show that Qf n 
is a Cauchy sequence in H“ p 0 (T) and the limit becomes a solution of <9 “m = A u+f. 


□ 


The following two results will be used later when we extend results proved for 
small T to the case when T is arbitrary. 

Lemma 5.2. Let u € H“ p 0 (T) and T < T. Then there exists u £ 0 (T) such 

that u(t) = u(t) for all t < T, and 

INIh“£(t) — (5.1) 

where Nq is from Lemma \5.1\ and is independent ofT. 

Proof. Denote / = d^u, and let u £ H“ p 0 (T) be the solution of 
uf = A u + (/ — Au)l t< f, t <T. 

Then by Lemma 15.11 

II^IIh“,'p(T) — N oUf- ^ u )^t<T IIl 9 ,p(T) < ^o||u|| H “.2(fy 

Next observe that for t <T, 

d?(u — u) = A u + f — Am — / = A(u — u). 

It follows from Lemma T5. II that u = u for t <T. The lemma is proved. □ 

Lemma 5.3. Let 0 < T < T and u,u £ H“ p0 (T). Assume that 

u(t) = u(t ) t <T (a.e.). (5-2) 

Then u(t) := u(T + t) — u(T + t) £ H“ p 0 (T — T). 

Proof. We take u ei,e2,e3 from the proof of Lemma |2~7TI (iii) which is defined as 

u ei ’ e2 ’ £3 (t,x)=r)(t)r) 3 (£ 3 x) / u{s,y)r] ltei (t-s)r) 2 ,e 2 {x-y)dsdy. 

J R d JO 

As shown before, for any e > 0 it holds that 

||m — u ei ’ £ 2 ' e 3 \\ K ; 2 {T) < £ 

if £i, £ 2 , and £3 are small enough. Hence we can take a sequence (a n , b n , c n ) so that 
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and 


Observe that 


||u — w an ’ 6 ”’ Cn || H j, 2 ( T ) —0 as n —> oo. 


u “n,b n ,c n ^ x j = Vi < (f + a„)AT 

due to ( 15 . 21 ) and the fact that 7yi i£l 6 C'“((£i,2ei)). Thus 

u„(i, x) := U a "-’ K ’ C " (f + t, x) - u a "’ fe "’ c " (f + t , x) 
is a defining sequence of u such that 

u n (0,x) = 0 and — u„(0,x) = 0 Vx £ R d , Vn £ N. 
Therefore the lemma is proved. 


□ 


Proof of Theorem 12.91 

Step 1. Denote / = b l u x i + cu + f(u). Then 

11/0) - /MIUp < n \W-v\\ h i + || f(u) - f(y)\\ Lp 

< e\\u - v 11 Hp + K\\u - u||l p + ||/(u) - /(u)||l p . 

Thus considering / in place of / we may assume 6 * = c = 0. 

Step 2. Let / = /o be independent of u. Assume that are independent of 
(t,x). In this case obviously we may assume a *- 7 = <5 U , and therefore the results 
follow from Lemma 15. II 

Step 3. Let f = fo- Suppose that Theorem 12.91 holds with some matrix a = 
x)) in place (a y (t, x)). We prove that there exists £o = £o(Ao) > 0 so that if 

sup |a y (f,x) — x)| < £q, 

(t,x) 

then Theorem l2.9l also holds for (a u (t,x)) with 2A/ in place of W- To prove this, 
due to the method of continuity, we only need to prove that ( 12 . 121 ) holds given that 
a solution u of (]2. 1 ID already exists. Note that u satisfies 

d?u = a l °u x i xj + /, f = f + ( a 13 - a lJ )u x i xj . 

Hence by the assumption, 

IMIh“’ p 2 (T) ^ -NolKy - 7 - a l °)u x i x j\\^ qp {T) + A r o||/||L 9 iP (T) 

< iVosup I a lJ - d l3 \ ||u|| H a,a (T) + ^oII/IIl^pCT)- 

Thus it is enough to take £ 0 = ( 2N 0 )~ 1 . 

Step 4. Let a 1 - 7 = a 1 - 7 (x) depend only on x and / = f 0 . In this case one can 
repeat the classical perturbation arguments to prove the claims. Below we give 
a detail for the sake of the completeness. As before we only need to prove that 
there exists a constant No independent of u such that ( 12 . 121 ) holds given that u is 
a solution. By Steps 2 and 3, there exists £o > 0 depending only on p , q, a, K, and 
T such that the theorem holds true if there exists any point xo £ R d such that 

sup |a y- (x) — a y- (xo)| < £o- (5.3) 
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Recall that a 13 is uniformly continuous. Let (io < 1 be a constant depending on e o 
such that 

\a lJ (x) - a lJ (y)\ < £ 0 /2, if \x-y\< 46 0 . 

Choose a partition of unity bn, n = 1, 2, • ■ ■ so that cj> n = b( a ~]o " ) for some x n £ 
R d and b £ 0)) satisfying 0 < b < 1 and b = 1 for \x\ < 1. Denote 

4>n — 4>{ X 2S^ )• Then <j) n = 1 on the support of </>„, and u n = <j) n u satisfies 

dn ~ dn (M , n)x i xi T fn 


where 

a l n = finCl^ix) + (1 - 4>n)a l:, {x n ), 
fn f4*n 2 a 3 U x i {(pn^ x j d ^ 'li{^(f> n ) x t x j . 

It is easy to check (a%) satishes (EH and (15.311 with x n in place of xq. By [ 19] 
Lemma 6.7] and Step 3, if p = q then for any t < T, 


I U IItut0.2 


9 


Y ||W>: 


n 11 ^ 0,2 




<N 0 Y II fn 


l L g,p(*) 


< N 


OO 

Y, [ll U a:( < / > ™)®llL g , p (t) + H' U (^«)®®llL,,p(t) + 


pW 




■mm 


pW 


(5.4) 


(5.5) 


The last inequality above is also from Lemma 6.7 of m- We emphasize that 
equivalence relation in ( 15 . 41 ) and inequality ( 15 . 51 ) hold in general only if p = q or 
only finite <p n are non-zero functions. Hence, if q ^ p then we take sufficiently 
large R,M > 0 so that X)«=i bn (x) = 1 on B R and vanishes for |a;| > 2 R, and 
the oscillation of a 13 on the complement of B R / 2 is less then eo/2. Denote bo = 
1 — Yln—i bn- Then one can repeat the above calculations and use the relation 
Hull*! < £r||u ||//2 + 7V||u||l to conclude that for all t < T 


ii-Hh-p) < ^n-iil,p( t ) + ^n/iil, P (t) 


< N 


< N 


f\[\s-r) 

J o IJo 

f(t-s) 

J o 


— l+a 


(Pt Q M(r )|| ff 2 + \\f(r)\\ L )dr 


ds + N Wf\\l gAt) 


■ 1+Q iMfe (s) ds 


' mf\\t qiP (T)> 


where Theorem 12.61 (iv) is used in the second inequality. Hence by a version of 
Gronwall’s lemma (see e.g. US] Corollary 2]) we get 

ll<-(T) ^ll/lll,p(T)- (5-6) 

From equation (12.111) . we easily conclude that 


Hbf^llL^pCT) < TV11 z/.|ljjo.2( r ) 


|Lq )P (T) 1 


and therefore (15.61) certainly leads to the a priori estimate. 

Step 5. Let /o = /, and a 13 be uniformly continuous in ( t,x ). As before, we 
only need to prove the a priori estimate ( 12 . 121 ) . 

Let No be the constant from Step 4 so that a priori estimate ( 12 . 121 ) holds when¬ 
ever a 13 are independent of t. Take £o from Step 3 corresponding to this Nq. We 
will apply Step 3 with a 13 = a l3 (to,x) for some to- 
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Take k > 0 so that 


la^(t,x) — (s,x)\ < eo/2 , if \t — s\ < 2 k. 

Also take an integer N so that T/N < re, and denote Ti = iT/N. By Steps 3 and 
4 applied with aN (t.x) = a l ' J (0, x). (12.121) holds with T\ and 2N$ in place of T and 
No respectively. Now we use use the induction. Suppose that the a priori estimate 
(12.121) holds for Tk < T with Nq independent of u. This constant N 0 may depend 
on k. Take u from Lemma 15.21 corresponding to T = Tk . Denote 

u(t, x) = (u — u)(T + t,x), f(t , x) = f(T + t , x). 


Then one can easily check that u satisfies 


d?u = a l ^(T + t, x)u x i x j + / + (a*- 7 '(T + t,x) — 5^)u x i x j (T + t,x), t <T — T. 

Due to Lemma IQ1 u is contained in H“’” 0 (T). Thus by the result of Steps 3 and 
4 with a u = a *- 7 (Tk ,x), we have 




/pm) 


< (2iV 0 ) 9 


fTk+i 


'T k 


\\f + (a ij - 6 ij )u x i x i\\ q L dt 


< 

< 


N\\u\\ 


+ Il/H£ 


Q 


P (T) < *11/11 


q 


(T)’ 


where the third inequality is due to m and the last inequality is from the as¬ 
sumption. Hence 




< U 


;p(r k+ 1 ) 






<*ll/r 




As the induction goes through, the a priori estimate (12.121) is proved. We emphasize 
that for each k the constant Nq varies, however for each Tk we use the result in 
Step 4 and therefore the choice of Tk+i (or the difference \Tk+i ~ Tfc|) does not 
depend on k, and therefore we can reach up to T by finite such steps. 

Step 6 . Let fo = f and = Ylk=i a k (*’ x )h T k-\Tk] (*)• If (■ = 1, (12.121) comes 
directly from Step 5. If l > 1, we use the induction argument used in Step 5. The 
only difference is that to estimate the solution on [T/., Tfc+i) we use the result of 
Step 5, in place of the result in Step 4. 

Step 7. The general non-linear case. We modify the proof of Theorem 5.1 of 
fl9] . For each u £ H“’p 0 (T) consider the equation 


d?v = a ij + /(«), t < T. 

By the above results, this equation has a unique solution v £ 0 (T). By denoting 

v = IZu we can define an operator 1Z : 0 (T) —> H ^ 0 (T). By the results for 

the linear case, for each t <T, 


< Ne q \\u - w ll^0;2( t ) + NKj\\u - v \\l Jq ' P (t) 
<N 0 ^\\u-v\\ q ^, ; +m f (t-s)- 1 +a \\u~v\\ q 2 {s) ds, 

J 0 
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where N\ depends also on e, and Theorem 12.61 (iv) is used in the last inequality. 
Next, we fix e so that 6 := No£ q < 1/4. Then repeating the above inequality and 
using the identity 

f {t- Sl y 1+a f { Sl - S2 y 1+a ■ ■ ■ f (s n _i-s„)- 1+ “ds n ---dsi 

Jo Jo Jo 

{T(a)}” 


T{na + 1) 


r 


we get 




m i e m ~ k (T a Ni) k 


^E 

k =0 

< 2 rn 0 rn 

— 2m^ 2 ll U “ (T)‘ 


|| u _ ^11 9 

T(fca + 1) 11 "K,p(T) 


(d~ 1 T a NiT(a)) 

max l - wr, - 7 \— 

k \ F(fca + 1) 


k \ 


"H ?; P 2 (T) 


For the second inequality above we use Y^k=o (jc'j = 2 m - It follows that if m 

is sufficiently large then 7 Z m is a contraction in H“p 0 (T), and this yields all the 
claims. The theorem is proved. □ 


6 . Kernels p and q 

6.1. The kernel p(t, x). In this subsection, we prove Lemma lrTT ih In other words, 
we introduce a kernel p(t, x) which is integrable with respect to x and satisfies 

Hp(t,m) = E a (-t a \t;\ 2 ). ( 6 .i) 

Let r( 2 :) denote the gamma function which can be defined (see [2] Section 1 . 1 ]) 
for z £ C \ {0,-1, —2,...} as 


T(z) = lirn 


\n. z 


n\n 


n->oo z(z + 1 ) • • • (z + n) 

Note that T(z) is a meromorphic function with simple poles at the nonpositive 
integers. From the definition, for z € C \ {0, —1, —2,...}, 

zT(z) = F(z + 1 ). ( 6 . 2 ) 

By (16.21) . one can easily check that for k = 0,1,2,..., 


Res 2 = _fcr(;s:) = lim (z + k)T(z) = 


(6.3) 


>-*• ' ' ' fc! ’ 

where Res 2 = _fcF( 2 ) denotes the residue of r(,z) at z = —k. It is also well-known 
(see. e.g. [2] Theorem 1.1.4]) that if 9?[z] > 0 and 3?[w] > 0 then 

r^rH 


/»oo /*1 

T(z)= f-'e-'dt, / t*- 1 (l-t)“- 1 dt = 

Jo Jo I\~ + ui) 

Using Stirling’s approximation (see [5] Corollary 1.4.3] or [161 (1.2.3)]) 
r(z) ~ V2^e^-^ )losz e- z , \z\ -> oo, 


(6.4) 


(6.5) 
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one can easily show that for fixed a G R, ([IBl (1.2.2)]) 

|F(a + i 6 )| ~ V2n\b\ a ~^e~ a ~^ | 6 | —> oo. ( 6 . 6 ) 


Let to, n, p, v be fixed integers satisfying 0 < m < p, 0 < n < v. Assume that the 
complex parameters ci,..., c v , di,..., d M and positive real parameters 71 ...., 7 „ . 
Si,... ,6/j, are given so that Pi D P 2 = 0 where 

p i ~ G C : j G {1, ■ • •, to}, k = 0,1,2,... 

p 2 := | 1 ~ ^ + fc G C : j G {1,.. ., n}, k = 0, 1 , 2, ... J . 

If either m = 0 or n = 0, then by the definition P\ fl P 2 = 0. For the above 
parameters, the Fox H-function H(r ) (r > 0) is defined as 


jjmn | (ci,7l) *** (*-n,7n) (^ro-f-1,7n-f-l) ***(c^,7i') 

^ (di, di) • * * (d m , d m ) (d m _(_i, d m _|_i) • * ■ (d^, d^). 

j_ r ut=i T ( d i +^ nE r(i ^ cj - jjz) _ 

27ri A rij= m +i r (! - dj - <s,-z) nj=n+i r ( c i+7j*) r 


(6.7) 


where the contour Z, is chosen appropriately depending on the parameters. Some 
special cases needed in our setting are specified below. In this article, we addition¬ 
ally assume that parameters ci,..., c„ and d \,..., d^ are real and 


E 'i/ E > °» E^~ E 'Vi + J2 6 i- E S > > °- 


( 6 . 8 ) 


3=1 


i=n -\-1 


. 7=1 


3=m+l 


Under (16.81) . we can choose the contour L of two different types. Hankel contour 
Lh is a loop starting at the point — 00 + ipi and ending at the point —00 + ip 2 
where p\ < 0 < P 2 , which encircles all the poles of Pi once in the positive direction 
but none of the poles of P 2 . Bromwich contour L v is a vertical contour, which go 
from 70 — ioo to 70 + ioo where 70 G R and leaves all the poles of Pi to the right 
and all poles of P 2 to the left. Braaksma [5] showed that the contour integral (16.71) 
makes sense along Lh and L v for r G (0,oo), and two integrals along Lh and L v 
coincide (see [161 Section 1.2]). Furthermore, if YTj =1 ~ J2i =1 7* > 0 then H(r) 
is an analytic function on (0,oo) and can be represented (see [T6J Theorem 1.2]) as 


H ( r ) = E E Res *= 


i=l k—0 


njll T ( d 3 + S j z ) UU r (! - C 3 - lj Z ) r _ z 

rij= m +i r(i — dj - Sjz) nj=„ + i r (cj+ Hz) r 

(6.9) 


where d^ := — (dj + k)/di G Pi are poles of the integrand in the contour integral, 
i G {1,..., to}, and k = 0,1, 2, • • •. Note that on the negative real half-axis the 
Mittag-Leffler function 


can be written as 


£«(*) = E 

fc =0 


T(aA: + 1) ’ 


a > 0 , 


( 0 , 1 ) 

(0,1) (0, a) 


E a {-x) = H“ 


(x > 0). 


( 6 . 10 ) 
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Indeed, by (16.911 and (16.31) 


H 


( 0 , 1 ) 

(0,1) (0, a) 


= £ 


(~i) fc r(i + fc) fc ^ (~x) k 

k\T(ak + 1) ' r(ak + 1) 


y: Res z= _fc 
fc=C 

OO 

= £ 


k =0 
oo 


r (z)r(i - z) 
r(l-a«) 

= E a {-x). 


k—0 ' ' fc=0 

We introduce some notion to introduce the asymptotic behavior of H{r). Let j s 
(1 < js < to) be the number such that 


p s := — = mm 

°3. 


(min 0 := oo) 


where the minimum are taken over all -y 2 - so that djo is a simple pole. Similarly let 
j c (1 < j c < to) be the number such that 


dj c 

Pc-= = mm 


L^J 


where the minimum are taken over all j 2 - so that djo is a pole with order n c > 2. 
Here n c denotes the smallest number of the orders for non-simple poles. 

The following result can be proved on the basis of (16.91) . See (HU Corollary 
1.12.1] for the proof. 

Theorem 6.1. (i) If p s < p c , then for r < 1 

\H(r)\ < Nr Ps ; 

(ii) if p s > p c , then for r < 1 

\H(r)\ < Nr pc \ lnrl” 0-1 . 

An upper bound of H(r) on [l,oo) is also well-known if n = 0 and m = p in 
(16.71) . See [HE Corollary 1.10.2] and [TSJ (2.2.2)]. 

Theorem 6.2. Suppose that n = 0 and m = p in (ED - Then for r > 1, 

\H(r)\ < 7Vr (A+1/2)w_1 exp {-ur]- 1/,JJ r 1/u }, (6.11) 

where 


A : ~£ d i £ Ci + 

l=i *=i 


v — p 


:= £ ^ - £ ^ r i- n v n a 

j=l i= 1 j=l i=l 


- 7 < 


Theorem 6.3. Suppose that n = 0 amd to = p in (16.71) . Then 
d 


dr 


— _^ _1 uri l + 1 0 
^H-l/2+l 


H(r ) = — 



(ci,7i) '• 

■ (0,1) 

7* 

(di,Si) 

(d^, S/f) (1,1) 


( 6 . 12 ) 


Now we define p(t, x) so that (16.11) holds. Fix a £ (0, 2) and let 


p(t,x):= 7T ^ 1^1 

1 




( 1 , a ) 

(|, 1 ) ( 1 , 1 ) 


= 7r-^H|- d 


r(s+ ^)r(i + 2) /i , 2 


27ri r(l + az) 


dz. 


(6.13) 
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One can easily check that (I 6 . 8 |l is satisfied. Hence we can take the Bromwich 
contour L = L v , which runs along from —7 — ioo to —7 + ioo, i.e. 

L v := {z € C : SR[ 2 ] = — 7 } 

where 

0 < 7 < min{l, —-—, —} if d > 2 and 0 < 7 <min{-,—} if d = 1 . 

4 a 2 a 

Note that 

{ p s = 1 if d > 3 

p c = 1 , n c = 2 if d = 2 

Ps = \ if d = 1. 

By Theorem 16.II and Theorem 16.21 

p(t, •) G Li(R d ), Vf>0. 

For the asymptotic behavior of the integrand in (16.1311 we use Stirling’s approxima¬ 
tion. Write z = —7 + ir, r G (— 00 , 00 ). Then by (16.61) . for p G (0, 00 ), t G (0,oo), 
and large |r|, 


r ( 2 + z )r(i + z) f 1 .—a. 2 

t p 


r(l + az) V 4 


1 


< N ( -t~ a p 


W 2 1 |r| Cl e _C2|r| , 


(6.14) 


where 

Therefore, 


. d — 1 7T . . 

ci := - 7(2 - a) + —, c 2 :=-( 2 -a). 


sup 

z(zL 


r(f + z)r(i + z) /1 _ 


T(l + az) V 4 


-r V < iv -mV 


(6.15) 


because for fixed a ^ {0, —1, —2, —3,...} the mapping b 1 —> T(a + i&) is a continuous 
function and does not vanish on R. Hence 


T(tj + z)T(l + z) /1 _ 


/ l T(l + az) V 4 


t “p 2 1 dz 


< 


r( 2 -7 + ^(1 -7 + ^)) (\.- a 2 

M p 


T (1 — a"f + i err) 


< Nt~ ai p 21 |V + ^ |r| Cl e- C 2 |r| dr| < VM a V 


dr 


27 


(6.16) 


We remark that (16.141) . (I6.15|) . and (16.161) hold for any 0 < 7 < min{l, ^}. 

Now we prove to . First assume d > 2. By the formula for the Fourier trans¬ 
form of a radial function (see (42) Theorem IV.3.3]), we have 


2^/2 /»oo 

= — d—r / P 

l£l 2 Wo 


-iH 2 ° 


-i“V 


(l,a) 

(f,l) ( 1 . 1 ) 


d4_i(|^|p)dp, 


where Jd _ 1 is the Bessel function of the first kind of order | — 1, i.e. for r G [0, 00 ), 

(_l) fc / r s2k-l+d/2 




k =0 


fc!r(fc + |) 


©' 
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It is well-known (e.g. [IS (2.6.3)]) that if m > —1 then 

Jm(t) = 

Thus 


0 (f m ), t —>• 0 + 
t —> oo. 


pOO pi p oo 

J o |p”^ + 27 Td_i(|^|p)|dp<y^ p^^dp + J^ p~^ +2l ~^dp < oo (6.17) 


since 0 < 7 < min{l 


d-1 j_ 
4 ’ a 


}. By the definition of the Fox H-function, 




p~2 H 


1 r 

27 ri ,/ 0 


f"V 


(l.«) 

(|, 1 ) ( 1 , 1 ) 


J4_i(|^|p)dp 


f r( 2 +4T (1 + *0 (l t -a p 2 ] dz 

I l r(l + az) V 4 


Jd_ 1 {\£\p)dp. (6.18) 


Combining (16.161) and (16.171) . we have 

_d r(4 +z)r(i + z) f 1 _ 


IQ JL 


r(i + az) 


t a p J*-i(\£\p) 


\dz\dp < oo. 


Thus we can apply Fubini’s theorem to (16.181) . Furthermore, since 

Vzgl, 


by using the formula QS (2-6.4)] 


pOO 

/ p--^ 2z Jd M\p)dp = 2 -i- 2 ^|l- 1 + 2 * 
Jo 


n-z) 

r (1 + zY 


we get 


r 

Jo 


p-^12 


f-y 2 


(i , a ) 

(|,i) ( 1 , 1 ) 


27ri 


P 2 2z Ji_x{\i\p)dp 


Jf_i(|£|p)dp 

F(| + ^)F (1 + z) (l ± _ c 


r(l + az) 
(t~ a \t;\- 2 )- z dz 


-t~ 


dz 


/£(i±7£(i+7£(-^)^_ altl _ 2w 

2 d /2 2tAJ l T(l + az)r(^+z) 

Kl *" 1 1 f r(l - z)T(z) 


2 d ! 2 27ri J_ L F (1 - az) 


C t a m~ z dz 


1 ^ 


2 d /2 


-H 


11 


*“17 


(o,i) 


l£l 




( 0 , 1 ) ( 0 , a) 

where the last equality is due to (16.101) . Therefore (16.11) is proved for d > 2. 
Next let d= 1. By (13. 10|) . 

„a -1 


l£| 2 ' 


pOO 

C [£ Q (-f“|£| 2 )])( S ) = / e~ st E a (-ne)dt = 

Jo 

Thus it suffices to prove 

r °° „q-i 

■)}](«) = J o e- st F{p(t,-)}(Z)dt=^-^. 


(6.19) 







































30 


ILDOO KIM, KYEONG-HUN KIM, AND SUNGBIN LIM 


By Theorems 16.II and 16.21 it holds that for each s > 0, 

/ / e~ st \p(t, x)\dxdt < oo. 

J 0 J R 

Hence by Fubini’s theorem, C [T { p(t , •)}] (s) = T{C [p(-, a;)] (s)}. Due to (16.411 and 

(EH. 


nOO 

£ [p(',z0] («) = / e~ st p{t,x)dt 
Jo 


= 7r 2 \x\ 


‘ /e- 4, f 

Jo 2tt1 Jl 


i f r(i + 2:)r(i + 2 ) /i _ 


27ri J L r(l + az) 


-t~ a x 


7r 2 HI 1 


r(i + ^)r(i + z) /1 


7 T 2 HI 1 S 1 


10 

— Z 


2 m J L T (1 + az) 

jj(J+ z )T(l + z)(Js a \x\ 2 ^ Z dz. 


e~ st t az dt 


dz 


dz 


dt 


2 ni 


Furthermore by |TBJ (2.9.19)], 

hj L n \ +z)ni+z) G s “ |x|2 ) iz 

'1 


= H; 


4« a |*| 2 


(|,1) (1,1) 


= 2 


/2U,x3/ 2 


^i/2(s a/2 |a;|), 


where K rj (z) is called the modified Bessel function of the second kinefl or the Mac¬ 
donald function and it satisfies (see [TJ 9.7.2]) 

37 T 


Hence we obtain 


K 1 / 2 (z) = ^z 1/2 e z , |arg~| < ^ 

£[?(■,*)] ( s ) = — ex p{—■ s “ /2 l a; l}, 


which obviously implies that 

T{C\p(;x)} («)}(£) = 


s a ./2 — l r 00 


e~ ixi exp{-s a/2 \x\}dx 


ga /2 — 1 2 s “/ 2 s a—1 


s“ + ici 2 s“ + iei 2 ' 

Therefore (16.191) is proved, and (16.11) holds. 

6.2. Representation of q(t,x) and K(t,x). Let a £ (0,2) and recall 


p(t,x):=w =|x| d Hi2 


- 4 t~lx\ 2 


( 1 , a ) 

(f,l) (1,1) 


By Theorems 16.21 and (16.121) . if t ^ 0 and x ^ 0 then p(t 7 x) is differentiable in t 
and lim t _ > .o+ p(t, x) = 0. Thus we can define 

q(t,x) := “ I , K(t,x) := JW). 

[Dt p(t,x), a €( 0 , 1 ) dt 


^It has also been called the modified Bessel function of the third kind. 
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In this subsection, we derive the following representations: 


q(t,x) =Tr-i\x\- d t a ~ 1 U 2 1 


-Tt~ a \x\ 2 


(a, a) 

(f,l) (1,1) 


and 


K(t,x) = n 2 |x| d t 1 H ^2 
We consider the Bromwich contour 


— t~ a I TI 2 

.4* (f,l) (1,1) 


d 1 


( 6 . 20 ) 


( 6 . 21 ) 


r(a +z)T(l + z) /1 2 

r(l + az) V4 11 


\dz\ ds 


L v := {z £ C : Sft[z] = - 7 }, 0 < 7 < min{l, -, -}. 

2 a 


Let /9 > 0. By (16.16II . 

iy^i 

<N\x \ 21 [ (t-s) / 3 - 1 (t-s)- a7 ds <00. 
Jo 

Thus by Fubini’s theorem and (10.41) . 

4 p(*, J (t ~ sf^pis, x)ds 

\t-sf~ 1 


v-i\x\- d rt 


m 


i rw + zni + z )n 


27 ri J L r(l + az) 


dz 


ds 


n-%\x\~ d r r(| + z)r(i + g) 
27 ri J L r(l + az) 


1 


r(/3) 


= 71-2 1x1-^ 


i f r(l + z)T(i + z) (1_ 


2niJ L T(l+P + az) \4 


( t-a) p - 1 s az da 
dz 


xr dz 


t~ a Ixl 2 


( 6 . 22 ) 


= 7 r-^|x|- d t / 3 H?° ±t-“Ixl 2 


(1 + /3, a) 

(|,1) (1,1) 


Therefore, (16.201) is proved for a £ (1, 2). 

Next we differentiate kernels with respect to t. Write z = —7 + ir. Following the 
method used to prove (16.141) and (16.151) . for any a £ (0, 2) and /3 £ [0, 2) we have 


r(f + z)r(l + z) /1 


T(l + /3 + az) V 4 


dt 


< 7 V(d, 7 ,/ 3 )t M-t Q |x| 2 ) [( 7 2 + I t | 2 ) 1 / 2 (1|t|>iI t | Ci P e C2,r| + l| T |<i) 


(6.23) 
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Hence the time derivative of the integrand in (16.2211 is integrable in z uniformly in 
a neighborhood of t > 0, and thus by the dominated convergence theorem, 


d tt20 
* Hl2 


\t~ a \x\ 2 


(1 + a, a) 

(|,1) (1,1) 


i d I f r(f + *)r(i + z) /1 _ 


2iridt 1 J L Y{l + ot + az) \4 

t^_ f T(%+z)T(l + z )__fl^_ 
27ri J L T(1 + a + az) 


-t~ a \x\ 2 


dz 


Using the relation 


az 


az [ -t~ a \x\ 2 


a 


dz. 


we obtain 

'1 


d xt20 

dt Hl2 

-Ll 

2iri 


t~ a \x I 2 


F(a + az + 1) r(a + az) T(a + az + 1) ’ 


(1 + a, a) 

(|, 1 ) ( 1 , 1 ) 


r(f + z)T(i + z) _ r(§ + *)r(i + z)\ f\ x _ 


r(a + az) 


r(a + az + 1) 


t a \x\ 2 ) dz 


= t~ l H?S 


Thus, 


\t~ a \x\ 2 


(a, a) 

(f,l) (1,1) 


- aH 


20 


it- a \x 1 2 


(1 + a, a) 

(|, 1 ) ( 1 , 1 ) 


D) a p{t,x) = —I?p{t,x) 


d_ 
cit 
d_ 
tit 


— 4 I I — d,aij20 

n 2 \x\ t H 12 


-t _ “|:r| 2 


= 7r -d / 2 |x| -d < at a_1 H 2 2 - t~ a \x 12 


1 j.a u tt 20 

^ ^ Hl2 


= 7r -d/ 2 | a; |-^«- 1 H 
Similarly, using the relation 


-f - “|d 2 


7 t~ a \x\ 2 


(1 + a , a) 

(f,l) (1,1) 

(1 + a, a) 

(f,l) (1,1) 

(1 + a, a) 

(|, 1 ) ( 1 , 1 ) 

(a, a) 

(|, 1 ) ( 1 , 1 ) 


az 


T(az +1) r(az)’ 


we have 


K ^)=di H - 


\t~ a \x I 2 


= 7r-S|x|-‘ , t- 1 H?§ 
Therefore (16.201) and (16.211) are proved. 


( 1 , a ) 

(i,o ( 1 , 1 ) 

Tt““d| 2 


(0,o) 

(1,1) (1,1) 
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6.3. Estimates of p(t,x ) and q(t,x). In this subsection we prove Lemma [331 h) 
and Lemma f3. 2 1 Since the case a = 1 is easier, we assume a ^ 1. 

By (16.711 and (16.91) with n = 0 and to = /i, 


H M ° 


r 


(ci,7i) 
(dr, 5 1) 


m 


= EE 


lim 

z->-dik 



{Cv, lu) 

{dfj,, dfj,) 
riik -1 / 


m oo 

= EE Ref w ifc 

i=l fc=0 


IIj=i r (dj + djz) 

n;=i r ( c j+7j*) 


(z - IIj=i r ( rf 3 + Sjz) _ \ ) 

(n ik - 1 )! n;=i r ( c t + T?^) r / J ’ 



where di k = —(di + k)/Si € P\ is a pole of the integrand in the contour integral 
and riik is its order for i = 1 ,..., p, and k = 0, 1 , 2, • • •. 

Let R := t _ “|a;| 2 , and denote 


ttUR) := H 


and 


1+fc+Z 2+fc+Z 


tt9 ( td\ _tt2+/c+Z 0 

ri k,l\ n ') •— n l+fe+Z 2+fc+Z 


-i? 


-R 


(1 , a ) 

(0,1) 

• (0,1) 

(1,1) 

(1,1) ■■ 

• (1,1) 

(a, a) 

(0,1) •• 

■ (0,1) 

(|,1) 

(1,1) •• 

• (1,1) 


where k, l = 0,1,2, ■ ■ •. Then by (16.121) . 


At 


= -x i \x\- d - 2 [dR p k0 (R) + 2B p kl (R)} 

< fV|x|- d_1 dHj |0 (J2) + 2Hj il (iZ) 


and 


D xl D xi (\x\- d Kl 0 (R)) 

= \D xi (-x>|- d - 2 [dHP 0 (R) + 2BP kil (R)]) | 

= {{-Sijlx]-^ 2 + (d + 2)x i x j \x\- d - i )[dB p k0 {R) +2R p kl (R)} 
+ 2x i xd\x\- d ~ 4 [dRl 1 (R) + 2R p k2 (R)}\ 

<N\x\- d - 2 J2\dKl l _ 1 (R) + 2Rl l (R)\. 

i=i 


Inductively, for any m = 3,4, • • • and k = 0,1, • ■ ■, we have 


D n 

^ X 


(|zrX,oCR))| <n\ 


E \ d Ki-^ R ) +m h w 
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Hence 


DT ( \x\~ d d^Ui° 2 


\ R 


(1 ) a ) 

(f,l) (1,1) 


\d?D?p(t,x)\ = 

k =1 

n m 

< N\x\- d ~ m t- n E | d Ki-d R ) + m U R ) 

k= 1 1=1 

Similarly, 

m 

(\x\- d Kl 0 (R)) \<N\x\- d ~ m J2\ d Ki-iW + m lA R ) 


I DZ 


and 


| d?D™q{t,x)\<Nt 


—n-\-a—1 


J2\ d 


k=1 


m ( i^i— 


,oM) 


< N\x\- a ~ m t- 


(6.24) 


r d “ ro ^" + ““ 1 EE + 2H M ( i? ) | ■ ( 6 - 25 ) 

k=Q 1=1 

Now we prove ()3.2[) and (|3.4|) . If i? > 1, by Theorem 16.21 

< NR^ +k+l)R2 ~ a) expj-cri?^} 

and 

H q k l {R) < NR^ +k+l+1 - a)/(2 ~ a) exp{-aR^}, 
where a = (2 — a)a a ^ 2 ~ a l Hence by (16.241) 

n m 

\d?D™p(t,x)\ < N\x\- d ~ m t- n EE R(?+k+t)/( 2-a) eX p{-CTi?^} 

fc=l /=0 

< N(t~ a ^ 2 ) d+rn t~ n exp{—(<r/2)t _2 =^ |^| 2 ^} 

»(d+m) 


< AT ^ n exp{—(<r/2)t 2 -“|a;| 2 -“}. 


Similarly, by (16.251) 


At D™q(t,x)\ < Nt 2 


(4+"0 


—n+a—1 


exp{-(o-/2)t 2 “°|x| 2 -“}. 


Thus (13.21) and (13.41) are proved. 

To prove (13.31) . we recall (16.21) . For k,l = 0,1,..., denote 




T(f + z){T(l + z)} k+l+1 T(f + z)T(l + z) 


,k+l 


T(1 + az){T(z)} k+l 


T(1 + az) 
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First we assume that d is an odd number. Then 0j( ; has simple poles at d\j = — 1 —j 
and d 2 j = — f — j for j = 0,1, 2, • • •. Due to (16.311 and (16.911 . for R < 1 we have 




2 oo 

*=i i=o 

2 oo 

= Inn ((« - d i:j )O p kl {z)R- 


i—l j =o 


= £(-w) 


j=o 


fc+ i (-i) 3 ' r (f-i-i) ^ 

j! T(1 — a — ja) 


00 " d (-i)j' r(i-f-j) 


3 = 0 


j! r(l - ¥ - aj) 




= (- 1 ) 


F(|-l) 


fc+l x v 2 


i?+ - 


fc+i 


r(i - f^) 


T(l-a) 

where O(R) is bounded in (0,1] by (16.51) and (16.21) . Hence 

\R p k 0 (R)\ < N(R + R 1 / 2 ■ l d=1 ). 

Moreover since 


——^-Ri + 0 (R){Rl +1 + R 2 ), 


d - 


fc+Z-1 


+ 2 - 


D 


fc+Z 


= 0 l = 1 , 2 ,..., 


(6.26) 

(6.27) 


it holds that for all l > 1 


\dH p k>l _ 1 (R) + 2& ktl (R)\<NR. (6.28) 

If d is an even number, Q p k z has a simple pole at dij = — 1 — j for 0 < j < | — 2 
and a pole of order 2 at d 2 j = — | — j for j = 0,1, 2, • • •. Hence by (16.31) and (16.91) . 
for R < 1 we have 


= E Res -« 


j=0 




1^712 T ^ ( Hes 2 —d 


-2 


I=o 

00 




E ,5g, (O - d «> e E ,,(»)«-■) + E Jw, ^ (O - <%) 2 @E,,(*)«-) 


3— 0 

4- 2 


J=0 


j=o 


fc+ * (-ip r(f-i-j) 


j! r(l-a(l+j)) 


4rfl 1+j • 1<M2 


™ d \ k+l (-i)i + 2 ^'- 1 

-L 


j'=o 


I?2+ilni? 


j! (3 + J “ I ) 1 r(l-^-ja) 


+ Res 2!=( i 
3=0 




R*+i 


= (-l) 


W r (l-1) R 1 , - 


d\ fc+i (-l)^" 1 In I? 


(|-i)! T(l-f) 


+ 0(R)(R2 + R 2 ), 
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where 0(R) is bounded in (0,1] due to (16.511 and (16.211 . Hence 

|H£ 0 (i?)|< (R + R\nR-l d = 2 ). (6.29) 

Moreover by (16.271) again, if l > 1 then 

\dRl l _ 1 (R) + 2Rl l (R)\<NR. (6.30) 

Therefore due to (16.2411 . (16.261) . (16.2811 . (16.291) . and (16.3011 . we obtain for any R < 1 
\d?p(t,x)\ < N\x\- d t~ n (R + R 1 ' 2 ■ l d=1 +RhiR- l d=2 ), 


and for any R < 1 and m £ N 

\d?D™p(t,x)\<N\x\- d ~ m t- n R, 


where N depends only on d, m, n, and a. Therefore Lemma 13.114 11 is proved. 
Next we prove (13.51) . For k, l = 0,1,..denote 


r(f+ ^){r(i + z)} fc + i + 1 r(f + 4T(i + *) k+l 
k ' l[Z) r(a + az){T(z)} fe +' r(a(l + z)) ~ ' 

We remark that 0^ t (z) does not have a pole at z = — 1 unless d = 2. First we 
assume that d is an odd number. Then 0| ; has a simple pole at dy = —2 — j and 
d 2 j = —| — j for j = 0,1, 2, • • •. By (16.311 and (16.911 . for R < 1 we have 


H K i? ) = EE Res - 

2—1 j =0 


2 oo 




E E E™ ((* ~ d b)*) 


2=1 j =0 


E(- 2 - i) 


/c+Z 


j=0 


(-l) 1+ j r(f-2-j) 
(1+i)! F(—a(l + j)) 


R 2+j 


i=o 


fc+i 


(_i)i r(i-f-j) 


= -(- 2 ) 


fc+i r(g-2) n2 

F(-a) 


j! r(o-^-jo) 

Zc+Z -p / -| d 


Rl+3 


- II- 


^\g Ri + Om R ^ + R% (6.31) 

1 o / 


where 0(i?) is bounded in (0,1] as before. Hence 

\Hl 0 (R)\<N(R 2 + R d/ 2 -l d=1}3 ). 

Furthermore, by (16.2711 and (16.311) if l > 1 then 

MH l ^R) + 2H^(J?)| < N(R 2 + R 3 ' 2 ■ l d=1 ). 


(6.32) 


(6.33) 


On the other hand, if d is even and greater than 2 then 0| ; has a simple pole 
at d\j = — 1 — j for 1 < j < | — 2 and a pole of order 2 at d 2 j = — f — j for 
j = 0,1,2, • • •. If d = 2, then 0^ ; has a simple pole at —1 and a pole of order 2 at 
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d, 2 j = —1— j for j = 1, 2, • • •. Thus by (16.3|1 and (16.91) . for R < 1 we have 


#-2 


KiW = E Res - 


=di 




ld>6 + Res z= _i 


i 


d =2 


7=1 

+ Res 2=c j 20 


'I 


' ld>4 + E! R es 3=d 2J 
7=1 

E ^. ((* - l d >6 + (-I) fc+i R 


®hR- 


• J-ci—2 


i=i 


7 OO 7 

+ ,5s. * (<* - « 2e h(*)- R ”) l>4 + E ** * (<*- d vfw>- 


7 = 1 


By the product rule of the differentiation, the above term equals 

-2 - ~rd 

l 2 

r(-aj) 

s/c+Z 


E(-1 -7)‘ +l • • r( l~ 1 V ) fil+i ' u » 

£-r j ! r(-aj) 


7=1 


+ (—!)” i? • ld =2 + Res z= d 20 i( z ) 


R 2 ■ ld>4 



ld>4 + E R es z=d 2 j 
7=1 

R.i+3 In R 

' r(a - f - aj) 


©k.l (*0 


i?^ 


T(# - 2) 


= -(-2) fc+ * • E? • ^>6 + (-!) fc+ ' R • ld=2 


+ Res z=l 


z=d 2 o 


T(-a) 

©fcd(~) 


R 2 • lrf>4 — —77 


cA fe+ ' (-1) 


-1 


V 


R 2 Ini? 

(f — 1)! r(a-f) 


L d>4 


0(i?)(i? 3 + i?s +1 (l + |lni?|)), 


where 0(i?) is again bounded in (0,1]. Hence 

|H« 0 (i?)| < N(R 2 + R ■ l d=2 + i? 2 Ini?. • l d=4 ) (6.34) 

and by (16.271) 

l dH fc,l- 1 (i?) + 2Hg | ,(i?)| < iV(i? 2 + i? 2 lni?- l d=2 ). (6.35) 

Combining (16.251) . (16.321) . (16.331) . (16.341) . and (16.351) . for any i? < 1 and nonnegative 
integer m we have 


\d?D?q{t,x)\ 

< N\x\~ d ~ rn t~ n+a ~ 1 

+ JV|o:|- d i- n+a - 1 


(i? 2 +i? 2 In i? • ld= 2 ) 

i? 1/2 • ld=l + R ■ ld =2 


T i? 2 In i? • ld =4 


1 


m=0} 


where iV depends only on d, m, n, and a. Therefore Lemma 13.21' i') is proved. 
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Finally we prove Lemma ET^ld i). that is for any t ^ 0 an x ^ 0, 

d?p = A p, ^ = A q- (6.36) 

To prove the first assertion above one may try to show that their Fourier transforms 
coincide. But due to the singularity of A p(t, •) near zero, we instead prove 




d?p(t, x)h{t)(j>{x)dxdt 



A p(t, x)h{t)cj)(x)dxdt 


(6.37) 


for any (j) G C£°(R d \ {0}) and h(t) G C™((0,T)). 

Note that t~ a \x\ 2 > c > 0 on the support of hcj), and thus (13.21) implies that 
^(t,x),dfp(t,x), and A p(t,x) are bounded on the support of hcj). Hence both 
sides of (16.371) make sense. By the integration by parts, 


[ dtp(t,x)hdt= f p(T — t,x)D^H(t)dt, 

Jo Jo 


where H(t ) := h(T — t). Recall that by Parseval’s identity 


[ fgdx = f (?f)(?g)dZ, V/, ff GL 2 (R d ). 

J R d Jn d 


Considering an approximation of / by functions in L 2 (R d ) one can easily prove 
that Parseval’s identity holds if / G Li(R <1 ) and g is in the Schwartz class. Hence 
the left side of (16.371) equals 



p(T — t, x)4>{x)dx 


d?H(t)dt 


0 l_JR d 


F; a (-(T-t)“|£| 2 )J-(</>)(£R 


d?H(t)dt. 


Observe that E((j))(^) —> 0 significantly fast as |£| oo and d^H{t) = 0 if t is 
sufficiently small. Hence by (13.11) we can apply Fubini’s theorem and show that the 
last term above is equal to 


f f E a {-(T-tT\tf)d?H(t)dt 
J R d JO 


nm)dt 


/R d Jo 

r T p 


E a (-t a mh(t)T(Am)dtd^ 


/ / p(t, x)h(t)A(j)(x)dtdx = / / Ap(t,x)h(t)(j)(x)dxdt. 

I o J n d Jo Jn d 


For the first equality above we use the integration by parts in time variable and 
the fact d^E a (—t a |f| 2 ) = — |£| 2 i? a (—f“|£| 2 ), and Parseval’s identity is used for the 
second equality. Therefore the first assertion of (16.371) is proved. 

Next we prove A q = ^. Note that due to (13.21) . if x ^ 0 then D x p(-,x) 
and D 2 p(-,x ) are bounded on (0,T) uniformly in a neighborhood of x. Hence, 
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if (0,1), 


A,=A si Mf 



— s)p(s, x)ds 

— s)p(s, x)ds 
s)Ap(s,x)ds = D 1 ~ a d?p 


Similarly, if a £ (1,2), 


dp 

dt' 


A q = A (J t “-V) = / t “ _1 (Ap) = /r 1 ^ 

Thus the second assertion of (16.361) is proved. 


dp 

dt' 
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